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Abstract. The two-dimensional Ising model with fixed magnetization is studied 
using Monte Carlo techniques. At the coexistence line, the macroscopic, extensive 
droplet of minority spins becomes thermally unstable by breaking up into microscopic 
clusters. Intriguing finite-size effects as well as singularities of thermal and cluster 
properties associated with the transition are discussed. 
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The two-dimensional Ising model has attracted much interest in the past. It is 
conceptually simple, and several of its non-trivial properties can be determined exactly 
However, despite numerous studies, the model seems to be not yet completely 
understood. 

In this Letter, we shall deal with the thermal stability of a droplet of minority spins 
in the two-dimensional Ising model with fixed magnetization, monitoring the transition 
from a compact single, macroscopic and extensive droplet at low temperatures to an 
ensemble of small clusters at high temperatures. Although related cluster equilibrium 
properties have been studied rather carefully in the usual, grand-canonical Ising model 
0, |J, this topic seems to have been largely overlooked. Only recently, the corresponding 
thermal behaviour of an adatom or vacancy island of monoatomic height on a crystal 
surface has been investigated Q. 

In particular, we consider a square lattice with L 2 sites and full periodic boundary 
conditions. Neighbouring spins, S^) = ±1, may interact ferromagnetically, with the 
coupling term —JSiSj, J > 0. We assume that N 2 , out of totally L 2 , spins are ' 
spins, with the magnetization 

M = 1 - 2(N 2 /L 2 ) (1) 

being conserved when varying the temperature, T. 

The resulting phase diagram in the (M, T)-plane is known to display a transition of 
first order between droplet and stripe phases || |6], 0. The transition takes place in the 
coexistence region, which is bounded by T (M ) describing the temperature dependence 
of the spontaneous magnetization of the standard Ising model. In the thermodynamic 
limit, Mq is determined by M 



l-(sinh2J/(fc B To))- 4 . (2) 



Another intriguing feature of the droplet phase in the coexistence region, T < To, 
will be discussed in the following, the thermal stability of the droplet of minority spins, 
i.e. '— ' spins for M > 0. 

At M > 1/2, the N 2 minority spins form a single square droplet in the ground 
state. When increasing the temperature at fixed magnetization, the largest cluster will, 
of course, shrink, but it may be, at sufficiently low temperatures, still extensive, with the 
number of cluster spins, thermally averaged, < N c >, being proportional to L 2 (oc N 2 ). 
Accordingly, the droplet size n c defined by 

n c =< N c > /N 2 } (3) 

will be non-zero in the thermodynamic limit, L, N — > oo, with M (N/L) being 
constant. 

However, at and above the coexistence line, T > T , the largest cluster is expected 
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Figure 1. Typical equilibrium configurations of the two-dimensional fsing model, at 
M = 0.68, below (a: k B T/J = 2) and above (b: k B T/J = 2.25) the cluster transition, 
simulating a lattice with 100 2 sites; '— ' spins are denoted by dark squares. 
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to be non-extensive @, see Fig. 1. Thence, one expects a 'cluster transition', T d at or 
below T . Indeed, assuming that, in the coexistence region, see Ref. 8, 

-M d = M S = M (4) 

where is the magnetization in the droplet, M s is the magnetization in the rest of 
the system, and M is given by Equation (2), one may easily obtain n c (T), at fixed 
magnetization M, as 

n c (T) = (1 + M )(M - M)/(2M (1 - M)) (5) 

Obviously, the droplet size will then vanish linearly in the reduced temperature t = 
|T -T|/T as T-^T @. 

To investigate the cluster transition and to check this relation, we did Monte Carlo 
simulations using a nonlocal spin-exchange algorithm [10 as well as an efficient cluster 
algorithm studying square lattices of various sizes, L 2 , with L ranging from 25 to 
950, at various magnetizations M, ranging from 0.68 to 0.98. 

As illustrated in Fig. 2, the droplet size n c (M, L, T) seems to approach Equation (5), 
as one considers larger and larger lattices, being presumably valid in the thermodynamic 
limit, L, N — > oo. This behaviour holds at all values of M, we studied. In addition, 
one observes rather interesting finite-size effects. 

First, the cluster transition is signalled, in a finite system, by the turning point of 
n c (T), at T S (L). For moderate system sizes, with L up to, say, 750, T S (L) is found 
to vary (almost) linearly in 1/L, for different values of M. A straightforward linear 
extrapolation would yield cluster transition temperatures close to, but somewhat below 
To. Accordingly, one may expect subtle finite-size corrections to the observed linear 
behaviour.- In this context, attention is drawn to the non-monotonic size dependence 
of n c (T) close to the coexistence temperature To, see Fig. 2. Obviously, rather large 
systems need to be considered in the vicinity of To to allow a reliable extrapolation to 
the thermodynamic limit. 

Another interesting finite-size effect is found by studying the size dependence of the 
largest cluster < N c > close to T . For small lattices, the droplet size scales like L x , with 
the (effective [|12|]) exponent x, at fixed temperature, being quite small, e.g. at M = 0.92 
x is seen to be around 0.5 to 0.7. Only for larger lattices, the cluster seems to become 
more compact, with x ~ 2. The crossover, at L c , is characterised by a pronounced peak 
in the effective exponent x, with the peak height and crossover length L c increasing as 
T approaches T . Obviously, this phenomenon is not predicted by Equation (5) which 
assumes a compact, extensive droplet, as seems to be correct in the thermodynamic 
limit. 

Previous evidence for the cluster transition is rather scarce. In the context of 
an adatom island on a crystal surface, the related transition was noticed 0, but it 
was not discussed in the framework of the phase diagram of the corresponding Ising 
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Figure 2. Droplet size n c as function of temperature for lattice sizes L ranging from 
125 to 750, at fixed magnetization M = 0.92. Error bars stem from averaging over 
at least N = 10 Monte Carlo runs with different random numbers. The dashed line 
corresponds to Equation (5). 



model. In that study, a non-trivial critical exponent describing the vanishing of the 
droplet size n c on approach to the cluster transition had been obtained, estimating 
the transition temperature from a straightforward linear extrapolation of simulational 
data for moderate system sizes, as discussed above.- A possible hint on the cluster 
transition might be hidden in a renormalization group study of Saito on quenches 
in Ising models with conserved magnetization, finding, however, a transition below the 
coexistence line, without mentioning of the thermal stability of the droplet.- Attention 
is also drawn to the work of Kertesz on droplet stability in Ising models in an external 
field [[14j] as well as related recent work on Ising cubes [15|, [16j and cluster shapes ||. 

It looks quite promising to analyse various cluster properties close to the cluster 
transition, in analogy to percolation studies ||. For instance, we find that the second 
moment of the cluster size distribution, excluding the largest cluster, shows a pronounced 
peak close to T s , indicating a divergence at the cluster transition in the thermodynamic 
limit. Furthermore, the specific heat also displays a pronounced maximum close to T s . 

In summary, we conclude that there seems to be, in the two-dimensional Ising model 
with conserved magnetization, a cluster transition at the coexistence line, at which the 
droplet of minority spins looses, in the thermodynamic limit, its extensivity. We think 
that the interesting cluster and thermal properties as well as finite-size effects close to 
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that transition deserve to be studied in more detail in the future, both in two and three 
dimensions. 
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